Three-body model calculations for 16 C nucleus 
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We apply a three-body model consisting of two valence neutrons and the core nucleus 14 C in 
order to investigate the ground state properties and the electronic quadrupole transition of the 16 C 
nucleus. The discretized continuum spectrum within a large box is taken into account by using 
a single-particle basis obtained from a Woods-Saxon potential. The calculated B(E2) value from 
the first 2 + state to the ground state shows good agreement with the observed data with the core 
polarization charge which reproduces the experimental B(E2) value for 15 C. We also show that the 
present calculation well accounts for the longitudinal momentum distribution of 15 C fragment from 
the breakup of 16 C nucleus. We point out that the dominant (ds/2) 2 configuration in the ground 
state of 16 C plays a crucial role for these agreement. 
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Nuclei far from the (3 stability line often reveal unique 
phenomena originating from a large asymmetry in the 
neutron and proton numbers. One of the typical exam- 
ples is a neutron collective mode, which is characterized 
only by the neutron excitation, with negligible contri- 
bution from the proton excitation. A recent calculation 
based on the continuum quasi-particle random phase ap- 
proximation (QRPA) has, in fact, predicted the existence 
of such neutron mode in the low-lying quadrupole exci- 
tation in 24 [1]. 

Recently, the electric quadrupole (E2) transition from 
the first 2 + state at 1.766 MeV to the ground state in 
16 C has been measured at RIKEN [2]. The observed 
B(E2) value (0.26±0.05 Weisskopf units) has turned out 
to be surprisingly small, as compared to the known sys- 
tematics in stable nuclei. On the other hand, a distorted 
wave Born approximation (DWBA) analysis for the 16 C 
+ 208 Pb inelastic scattering indicates a large enhance- 
ment of the ratio of the neutron to proton transition am- 
plitudes, M n /M p =7.6±1.7 [3], that is considerably larger 
than the isoscalar value, N / Z—l.&I . A similar value for 
M n /M p was also obtained from the inelastic proton scat- 
tering from the 16 C nucleus [4]. These experimental data 
suggest that the first 2+ state in 16 C is a good candidate 
for the neutron excitation mode. 

There have already been various theoretical calcula- 
tions for the structure of 16 C nucleus [5-10]. Except for 
a recent microscopic shell-model calculation [10], how- 
ever, they all fail to reproduce the anomalously hindered 
E2 transition. For instance, Suzuki and his collaborators 
have solved a n+n+ 14 C three-body model and found that 
the E2 strength is overestimated by a factor of about 2 
if the same core polarization charge is employed as that 
used to describe the 15 C nucleus. A similar overesti- 
mation of B(E2) value was found also in the antisym- 
metrized molecular dynamics (AMD) calculation [5] as 
well as in the deformed Skyrme Hartree-Fock calculation 
[6]- 

In this paper, we apply a three-body model with a 
finite-range n-n interaction [8, 9] to describe the ground 
and excited states in the 16 C nucleus. We employ the 



single-particle (s.p.) basis obtained from a n- 14 C Woods- 
Saxon potential to diagonalize the three-body Hamilto- 
nian. The continuum s.p. spectrum is discretized in a 
large box. Notice that the effect of continuum couplings 
can be properly accounted for with such s.p. basis [11]. A 
similar three-body model with a density-dependent con- 
tact interaction has successfully been applied to describe 
the structure of Borromean nuclei [12-16]. In Refs. [8, 9], 
Suzuki et al. adopted the correlated Gaussian basis to di- 
agonalize a similar three-body Hamiltonian for 16 C. How- 
ever, it remains an open question whether the correlated 
Gaussian basis is efficient enough to take into account 
the continuum couplings. Therefore, our study can be 
considered as a complement to the previous studies in 
Refs. [8, 9]. 

Assuming that the effect of core excitation on the low- 
lying spectrum of the 16 C nucleus is negligible [8, 17], we 
consider the following three-body Hamiltonian: 



H = h(l) + h(2) + V n . 



Pi P2 



(1) 



where m and A c are the nucleon mass and the mass num- 
ber of the inert core nucleus, respectively, h is the s.p. 
Hamiltonian for a valence neutron interacting with the 
core. The diagonal component of the recoil kinetic en- 
ergy of the core nucleus is included in h, whereas the 
off-diagonal part is taken into account in the last term in 
the Hamiltonian (1). We use a Woods-Saxon potential 
for the interaction in h, 



exp 
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V nC (r) = (v + Vi s (l-s)--p\ \l 
\ r ar j 
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where R = tqA c . The parameter sets for the Woods- 
Saxon potential which we employ in this paper are listed 
in Table I. The sets A, B, and C were used in Ref. [8, 9], 
while the set D was used in Ref. [17] in order to discuss 
the role of particle- vibration coupling in the 15 C nucleus. 
These parameter sets yield almost the same value for the 
energy of the 2si/ 2 state, €2s 1/2 ~ —1.21 MeV, and of the 
the ld 5 / 2 state, e lds/2 ~ -0.47 MeV. 
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TABLE I: Parameters for the Woods-Saxon neutron-core po- 
tential, V n c, in Eq. (2). 



Set 


V (MeV) V ls 


(MeV fm 2 ) 


r (fm) 


a (fm) 


A 


-50.31 


16.64 


1.25 


0.65 


B 


-50.31 {1=0) 


31.25 


1.25 


0.65 




-47.18 (I / 0) 








C 


-51.71 


26.24 


1.20 


0.73 


D 


-44.41 


31.52 


1.27 


0.90 



In our previous works [15, 16], we used the density- 
dependent delta force [12, 13] for the interaction between 
the valence neutrons, V nn . However, we here use the 
same finite-range force as in Ref. [8] in order to compare 
our results with those of Refs. [8, 9]. That is the singlet- 
even part of the Minnesota potential [18], 

V nn (r u r 2 ) = v a e-Mn-r,)' + ^ e - bl{ri -r 2 )^ (3) 

with v =200 MeV, 6 =l-487 fm" 2 , and 6 1= 0.465 fm" 2 . 
Following Ref. [8], we adjust the value of vi for each 
parameter set of the Woods-Saxon potential so that the 
ground state energy of 16 C, E gs — —5.47 MeV, is repro- 
duced. 

TABLE II: Properties of the ground and the second + states 
obtained with several parameter sets for the neutron-core po- 
tential, V„c- Pas and P dd are the probability for the [(2s!/ 2 ) 2 ] 
and [(W5/2) 2 ] components in the wave function, respectively. 
Ps=o is the probability of the spin-singlet (S = 0) compo- 
nent in the ground state. E + is the excitation energy for 

the second + state in the unit of MeV, while r( 16 C) is the 
root-mean-square radius of the 16 C nucleus in the unit of fm. 
The experimental values are E Q + (exp) =3.00 MeV and r( 16 C; 
cxp)=2.64 ± 0.05 fm [19], respectively. 



Set 


Pss(g-S 


) Pdd(g-S 


) Ps=o 


r( 16 C) 




Pss(0t) Pdd(0t) 


A 


0.184 


0.699 


0.784 


2.56 


2.32 


0.755 


0.201 


B 


0.177 


0.711 


0.746 


2.56 


2.35 


0.775 


0.187 


C 


0.183 


0.696 


0.768 


2.57 


2.39 


0.763 


0.196 


D 


0.206 


0.633 


0.808 


2.64 


2.48 


0.733 


0.221 



The three-body Hamiltonian (1) is diagonalized by ex- 
panding the two-particle wave function 5 , (ri,r 2 ) with 
the eigenfunction <j) n ij m of the s.p. Hamiltonian h, where 
n is the radial quantum number. The continuum s.p. 
states are discretized with a box size of -ffbox = 30 fm. 
We include the s.p. angular momentum l\ and Z 2 up to 5, 
and truncate the model space of the two-particle states 
at ei + e 2 = 30 MeV, where e is the s.p. energy of the 
valence particle. We have checked that the results do not 
significantly change even if we truncate the model space 
at 60 or 80 MeV, as long as v\ in Eq. (3) is adjusted for 



each model space. In the diagonalization, we explicitly 
exclude the lsi/ 2 , lp3/ 2 , and lpi/ 2 states, which are occu- 
pied by the core nucleus. The results for the ground state 
and the second + state are summarized in Table II. The 
parameter set dependence is small, although the set D 
reproduces the excitation energy of the second + state, 
E Q +, and the root-mean-square (rms) radius of the 16 C 

nucleus, r( 16 C), slightly better than the other parameter 
sets. The latter quantity is calculated as [12-14] 

(4) 

where A = (ri + r 2 )/2 and p = r\ — r 2 . Following Refs. 
[8, 9], we take 2.35 fm for the rms radius of the core 
nucleus, Vv^T- We find that the rms radius of 16 C is 
well reproduced in the present calculations. 

We notice that our results are considerably different 
from those of Refs. [8, 9] concerning the probability for 
the [(2s!/ 2 ) 2 ] and [(ld 5 / 2 ) 2 ] components in the wave func- 
tion, which are denoted by P ss and Pdd in Table II, re- 
spectively Our results show that the ground state of 16 C 
mainly consists of the [(lds/ 2 ) 2 ] configuration, while the 
second + state is dominated by the [(2s!/ 2 ) 2 ] configu- 
ration. This is in contrast to the results of Refs. [8, 9], 
which show the dominance of the [(2s!/ 2 ) 2 ] component 
in the ground state. As a consequence, we also obtain a 
smaller value of the spin-singlet probability, Ps=o, than 
in Ref. [9]. Notice that the d-wave dominance was sug- 
gested from the analyses of longitudinal momentum dis- 
tribution for the one- neutron knockout reaction of 16 C 
[20, 21]. We will discuss the longitudinal momentum dis- 
tribution later in this paper. 

It is worthwhile to consider a simple two-level pair- 
ing model consisting of the 2s!/ 2 and lc?5/ 2 s.p. lev- 
els in order to illustrate how the [(ld 5 / 2 ) 2 ] configuration 
becomes dominant in the ground state of 16 C. If there 
were no interaction between the valence neutrons, the 
ground state wave function would be the pure [(2si/ 2 ) 2 ] 
state, since the s.p. energy for the 2si/ 2 state is lower 
than that for the ld 5 / 2 state (e 2si/2 = —1.21 MeV and 
eid 5/2 = —0.47 MeV) . If one assumes a delta interaction, 
V nn = —g5(ri — r 2 ) between the valence neutrons, the 
diagonal matrix element of the Hamiltonian reads [22] 

H u = 2t t - g 2j Q + 1 hi , (5) 

where In is the radial integral for the configuration 
(i) 2 . Therefore, the pairing interaction influences the 
[(ld 5 / 2 ) 2 ] configuration more strongly than the [(2s!/ 2 ) 2 ] 
configuration by a factor of 3 when the radial integral 
In is similar to each other. If we choose the strength g 
so that the ground state energy is reproduced within the 
two-level model, we find <?=1005 MeV-fm~ 3 for the pa- 
rameter set D. This leads to the diagonal matrix element 
of H lt = -3.70 MeV for i = 2s 1/2 and H H = -4.85 MeV 
for i = ld 5 / 2 , lowering the [(l<i5/ 2 ) 2 ] configuration in 
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energy. Taking into account the off-diagonal matrix ele- 
ment and diagonalizing the 2x2 matrix, wc find P ss =0.26 
and Pdd=0.74 for the ground state, which are very close 
to the results shown in Table II. 

The upper panels of Figs. 1 and 2 show the two- 
particle densities p 2 [r\ ,r 2 ,0) for the ground state and 
the second + state, respectively. These are obtained 
with the Minnesota potential and the parameter set D 
for the s.p. potential. In order to facilitate the pre- 
sentation, we set n = r 2 = r and multiply the weight 
factor of 87r 2 r 4 sin# [15]. Despite that P ss and Pdd are 
considerably different, we obtain similar density distribu- 
tions to those in Ref. [9]. In particular, we observe sim- 
ilar "di-neutron" and "cigar-like" configurations in the 
ground state, as well as "boomerang" configuration [9] 
in the second + state as in Ref. [9]. The lower panels 
of Figs. 1 and 2 show the angular densities p(9) ob- 
tained by integrating the radial coordinates in the two- 
particle density [15]. It is multiplied by a weight factor of 
27rsin0. As a comparison, we also show the angular den- 
sities for the pure [(2si/ 2 ) 2 ] and [(ld 5 / 2 ) 2 } configurations 
by the dotted and dashed lines, respectively. They are 
given by p{9) = l/An for the [(2s!/ 2 ) 2 ] configuration and 
p{9) = 3/4tt -(5/4- cos 4 9 + 3/20) for the [(ld 5/2 ) 2 ] config- 
uration. As we see in the figures, the angular density for 
the ground state is close to that for the pure [(ld 5 / 2 ) 2 ] 
configuration, while the angular density for the second 
+ state is close to that for the pure [(2si/ 2 ) 2 ] configura- 
tion, being consistent with the calculated values for Pdd 
and P ss listed in Table II. 

We next discuss the quadrupole excitation in 16 C. Ta- 
ble III summarizes the results of the present three-body 
model for the first 2+ state. The energy of the 2 + state 
is well reproduced with this model, especially with the 
parameter set D. As compared to the results of Refs. 
[8, 9], the probabilities for the [2si/ 2 lci5/ 2 ] and [(ld 5 / 2 ) 2 ] 
components, denoted as P s d and Pdd, respectively, are 
comparable to each other in our calculation, whereas P s d 
is much larger than Pdd in Ref. [9] . 

In order to calculate the E2 transition strength, we 
introduce the core polarization charge, e po i. The E2 op- 
erator Q 2fl in the present three-body model then reads 
(for p=0) [8, 9], 



Z r A 2 -2A + 2 » ^ 2 

e H - e pol J 2^ r i l^oK) 



A 2 



A 2 



5 ( Z c n A-l \ .„ 

^I^e- 2 A2 e pol I (2ziz 2 - x x x 2 - 2/i2/2j,(6) 

where A — A c + 2. The value of the core polarization 
charge which is required to fit the experimental B(E2) 
value in the 15 C nucleus, B(E2; 5/2+ -» 1/2+) = 0.97 ± 
0.02 e 2 fm 4 [23], is listed as cp | in the fifth column in 
Table III. Notice that these are significantly smaller than 
that obtained with the harmonic vibration model of Bohr 
and Mottelson, esM = 0.55 for 16 C, which, however, does 
not include the effect of loosely-bound wave functions 
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FIG. 1: (Color online) (the upper panel) The two-particle 
density for the ground state of 16 C obtained with the param- 
eter set D for the single-particle potential as a function of 
ri = ?"2 = r and the angle between the valence neutrons, 8. It 
is weighted with a factor of 8n 2 r 4 sin 8. (the lower panel) The 
corresponding angular density weighted with a factor 2ir sin 8. 
The solid line is the result of the three-body model calcula- 
tion with the Minnesota potential. The dotted and the dashed 
lines are for the pure [(2s 1/ / 2 ) 2 ] and [(lds/2) 2 ] configurations, 



respectively. 



(see Eq. (6-386b) in Ref. [24]). For a loosely-bound 
state, the polarization charge may be modified as 



3i? 2 /5 



{l2j2\r 2 \hji 



(<) 



where R = 1.2A 1 / 3 fm and (hj2\r 2 \hji} is the radial ma- 
trix element between the s.p. states (hji) and (l 2 j 2 ) (See 
Eq. (6-387) in Ref. [24]). With the set D, we obtain the 
ratio 1 (lds/ 2 \r 2 |2s!/ 2 ) =0.205 as a reduction factor for 
the polarization charge for the transition from ldg/2 to 
2s\j 2 s.p. states. This leads to e po i=0.113, which is con- 
sistent with e pol shown in Table I. A similar small value of 
polarization charge has been obtained also with the self- 
consistent Hartree-Fock (HF) + particle-vibration model 
[25]. The calculated B(E2) value for 16 C with e pol is 
listed in the sixth column in Table III. In contrast to the 
previous calculations with the three-body model [8, 9], 
which overestimated the B(E2) value for 16 C with e pol , 
our calculations reproduce well the experimental B(E2) 
value. We notice that the small value of P ss and P s d in 
our wave functions is responsible for good agreement with 
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FIG. 2: (Color online) Same as Fig. 1, but for the second + 
state of 16 C. 



TABLE III: Properties of the first 2+ state of 16 C obtained 
with several parameter sets for the neutron-core potential, 
Vnc- Psd and Pdd are the probability for the [(2s!/ 2 ) x (W5/2)] 
and [(Ids/2) 2 ] components in the wave function, respectively. 
E 2 + is the excitation energy in the unit of MeV, while B(E2) is 
the electric quadrupole transition strength from the 2 + state 
to the ground state, in the unit of e 2 fm 4 . The experimen- 
tal values are E 2+ (exp) = 1.77 MeV and B(E2; exp)= 0.63 
± 0.27 e 2 fm 4 [2], respectively. ep 0l is the core polarization 
charge which reproduces the experimental B(E2) value for 
the 15 C nucleus within the n- 14 C model, whereas e^ ol takes 
into account the mass number dependence according to Eq. 
(11). 



Set 


E 2+ 


Pad 


Pdd 


^pol 


B(E2; 4 ol ) 


11 

^pol 


B(E2; e^i) 


A 


1.26 


0.392 


0.504 


0.162 


0.972 


0.145 


0.808 


B 


1.33 


0.400 


0.515 


0.160 


0.937 


0.144 


0.781 


C 


1.34 


0.402 


0.500 


0.153 


0.956 


0.137 


0.797 


D 


1.63 


0.472 


0.406 


0.122 


1.074 


0.109 


0.899 



the experimental B(E2) value. For the parameter set D, 
the E2 matrix elements between various two-particle con- 
figurations are estimated to be, 



(2s 1/2 ld 5/2 



IQ2o|(25i/ 2 ) 2 )- -1.087 efm 2 , (8) 



(2s 1/2 14/2IQ20I [(ld 5 /2 



,2 1,7=0 



[(Ids/: 



\2l J=2\ 



>|[(l<*5/5 



\2l J=0 



) = -0.627 efm 2 , (9) 
) = -0.811efm 2 10) 



Thus, the largest matrix element is the one between 
the (2S1/2) 2 configuration in the ground state and the 
[2S]/2 1^5/2] configuration in the 2 + state, although the 
other two matrix elements have substantial contribu- 
tions. Naturally, a small value of P ss and P s d leads to 
a small B(E2) value, which is desired in order to repro- 
duce the experimental data. A further improvement of 
the calculated value of B(E2) can be achieved if the mass 
number dependence of polarization charge is taken into 
account. In Ref. [26], the result of the HF+particle- 
vibration model for the core polarization charge of car- 
bon isotopes was parameterized as, 



e pol = 0.82- 



0.25- 



N - Z 



0.12 



Z N -Z 

0.36- — 

A A 



(11) 

This formula leads to the ratio of e po i for 16 C to that for 
15 C to be e po i( 16 C)/e po i( 15 C) =0.897. The polarization 
charge which is scaled by this factor from e pol is denoted 
by e^ ol in Table III. We find that the calculated B(E2) 
values with ep 0l agree remarkably well with the experi- 
mental value within the experimental uncertainty. 

We now discuss the longitudinal momentum distribu- 
tion of 15 C fragment in the breakup reaction of 16 C nu- 
cleus. For this purpose, we calculate the the stripping 
cross section in the eikonal approximation [27-30] . That 
is 



29, 30], 

da- n _ 

dpz 



1 



2nh 21 



—T 
+ 1 z - 



d 2 b n [1- |5„(6 n )| 2 



d 2 r 1 _\S c (b c )\ 



dze~^ z l h g l3 {r)Y lm (r) 



(^2) 



where gij(r) is the radial part of the spectroscopic am- 
plitude given by <1^ m ( 15 C)|* 3S ( 16 C)} = 5j; (r-)^-™(r) 
with 3^z_ m (r) being the spinor spherical harmonics. b n 
and b c are the impact parameters for the neutron and the 
core nucleus, respectively. They are related to the rela- 
tive coordinate between the neutron and the core nucleus, 
r = (r±,z), by b n = b c + r±. In the eikonal approxi- 
mation, the S-matrix is calculated as S(b) = cxp(2ix(fo)) 
with 



dz V(b, z), 



(13) 



where v is the incident velocity and V(b, z) is an optical 
potential between a fragment and the target nucleus. 

Figure 3 compares the eikonal approximation for the 
breakup reaction 16 C + 12 C ^ 15 C + X at E = 83 
MeV/nucleon with the experimental data [21]. We use 
the optical potential of Comfort and Karp [31] for the 
neutron- 12 C potential. The optical potential between the 
15 C fragment and the target is constructed with the sin- 
gle folding procedure using the 14 C density given in Ref. 
[32] and a s.p. wave function for the valence neutron 
for a specified final state of the fragment nucleus. As is 
often done, we assume that the cross sections for diffrac- 
tive breakup (i.e., elastic breakup) behave exactly the 
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FIG. 3: The longitudinal momentum distribution of the C 
fragment from the breakup reaction of 16 C on 12 C target at 83 
MeV/nucleon, obtained with the parameter set D. The contri- 
bution from the l/2 + and 5/2 + states of C are shown by the 
dotted and the dashed lines, respectively. The experimental 
data are taken from Ref. [21]. 



same as the stripping cross sections as a function of lon- 
gitudinal momentum, and thus we scale the calculated 
cross section (12) to match with the peak of the exper- 
imental data. The dotted and dashed lines in Fig. 3 
are contribution from the lsi/2 and 2d 5 / 2 states of the 
fragment C nucleus, respectively. These are added in- 
coherently to obtain the total one neutron removal cross 
section, which is denoted by the solid line. Our result re- 
produces remarkably well the experimental longitudinal 
momentum distribution of the 15 C fragment in the range 



of -200MeV/c <p z < 200MeV/c. 

In summary, we have applied the n-n- C three-body 
model in order to investigate the properties of the 16 C 
nucleus. We diagonalized the three-body Hamiltonian 
with the finite range Minnesota potential for the interac- 
tion between the valence neutrons. As the basis states, 
we adopted the single-particle states obtained from the 
Woods-Saxon potentials, in which the continuum spec- 
trum is discretized within the large box. With this model, 
the experimental data for the root-mean-square radius, 
the B(E2) value from the first 2+ state to the ground 
state, and the longitudinal momentum distribution of the 
15 C fragment from 16 C breakup are all reproduced well. 
In particular, we have succeeded to reproduce the B(E2) 
value for both 15 C and 16 C nuclei simultaneously using 
the core polarization charge which is consistent with the 
one obtained with the particle- vibration coupling models. 
The calculated probability of the (lds/2) 2 configuration 
in the ground state wave function of 16 C is about 70% 
while that of the (2s!/ 2 ) 2 configuration is about 18%. 
These values are close to those extracted from the anal- 
yses of the experimental longitudinal momentum distri- 
bution. 
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